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SECTION  I 
INTRODUCTION 


It  is  now  generally  accepted  that  flight-by-flight  fatigue 
testing  is  more  representative  of  service  experience  than  block 
type  testing.  Further,  It  appears  that  load  randomization  within 
a  flight  corresponds  reasonably  well  to  flight  time  history 
measurements.  The  purpose  of  this  report  Is  to  provide  a  computer 
technique  to  transition  from  exceedance  data  for  various  phases  of 
flight  to  a  randomized  fllght-by-f light  loading  sequence. 

The  method  is  based  on  the  assumption  that  the  exceedance 
function  for  stress  at  a  control  point  of  the  aircraft  has  been 
stepped  (  i.e.  defined  in  terms  of  the  number  of  occurrences  at  a 
finite  number  of  stresses).  One  approach  for  obtaining  the  stepped 
exceedance  function  is  described  in  reference  (1).  These  data  are 
derived  for  take-off  taxi,  flight  (including  all  of  the  mission 
interval  divisons  such  as  ascent,  cruise,  combat,  loiter,  and  descent), 
landing  impact,  touch-and-go  landings,  and  landing  taxi  and  arranged 
for  input  to  the  computer  routine.  In  addition  temperature  exceedance 
functions  may  be  placed  in  the  program  to  provide  a  temperature  cor¬ 
responding  to  each  stress  occurrence. 

It  is  assumed  that  the  stresses  and  temperatures  are  randomized 
within  each  phase  of  flight.  The  ordering  of  missions  flown  within 
a  lifetime  may  be  predetermined  or  determined  randomly.  In  addition, 
some  freedom  Is  allowed  in  selecting  the  flight  time  in  a  mission 
randomly.  This  feature  is  useful  in  the  case  where  the  exceedance 
data  are  derived  from  flight  recorder  data. 

In  all  flight  phases,  the  random  sampling  is  done  without 
replacement.  As  a  consequence,  the  exceedance  function  derived  from 
the  fllght-by-f light  load  sequence  is  identical  with  the  Input  ex¬ 
ceedance  function. 

The  computer  program  along  with  a  sample  problem  is  included. 

This  program  terminates  at  the  point  of  writing  the  fllght-by-f light 
sequence  of  loads.  To  be  of  practical  value,  the  output  of  this  routine 
will  have  to  put  on  cards  or  tape  (whatever  is  required  by  the  specific 
requirements  of  the  test  set-up). 
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SECTION  II 


DESCRIPTION  OF  METHOD 


Based  on  flight  recorder  data,  or  for  a  new  design,  the  mission 
profiles,  the  stress  spectrum  for  positive  load  factors  and  the  stress 
spectrum  for  negative  load  factors  can  be  computed.  The  method  for 
doing  this  was  developed  in  Reference  1.  These  results  when  combined 
with  gust,  landing,  and  taxi  loads  can  be  used  to  generate  the  fatigue 
test  loading  sequence. 

There  are  two  possibilities  that  will  be  considered  for  this 
sequence.  These  are 

Case  1.  -  Stress  spectra  derived  from  flight  recorder  data 

For  this  case  the  basic  assumptions  are 

a.  Loading  randomized  within  a  specified  weight 
interval  and  loadings  per  unit  time  fixed  within 
an  interval 

b.  Weight  Intervals  treated  as  nonrandom  within  a 
given  mission  type 

c.  Mission  type  treated  as  random  or  nonrandom  within 
aircraft  life 

d.  Hours  per  flight  are  randomized 

e.  Temperature  treated  as  random  within  a  given 
weight  interval 

Case  2.  -  Stress  spectra  derived  from  the  mission  profiles 

The  basic  assumptions  are 

a.  Loading  randomized  within  a  mission  interval  and 
loadings  per  unit  time  fixed  within  a  interval 

b.  Mission  Intervals  treated  as  nonrandom  within  a 
given  mission  type 

c.  Mission  type  treated  as  random  or  nonrandom  within 
aircraft  life 

d.  Hours  per  flight  determined  from  the  appropriate 
mission  profile  or  randomized 

e.  Temperature  treated  as  random  within  a  given 
mission  interval 
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For  both  of  these  cases  the  following  assumptions 
will  be  made 


(a)  The  sampling  of  the  stress  spectra  will  be 
done  without  replacement 

(b)  When  treated  as  a  random  process  the  mission 
type  will  be  sampled  without  replacement 

The  flight-by-flight  fatigue  test  loading  sequence 

for  a  specified  mission  is  developed  in  the  following 

order: 

(a)  Aircraft  in  take-off  configuration  resting  on 
the  landing  gear 

(b)  Take-off  taxi  load  cycles  applied  for  the 
time  appropriate  to  this  phase.  It  is 
assumed  that  the  first  load  is  a  maximum, 
then  a  minimum  followed  by  a  maximum  etc. 
Further,  It  is  assumed  the  last  load  is  a 
minimum. 

(c)  The  flight  load  cycles  are  applied  mission  or 
weight  (as  appropriate)  interval  by  Interval  for 
the  entire  flight.  The  first  loading  in  flight 
is  a  maximum.  This  is  followed  by  a  minimum 
(either  from  the  population  of  one  g  load  or 
the  population  of  less  than  one  g  loads).  The 
maximum-minimum-maximum  etc.  sequence  is  con¬ 
tinued  throughout  the  segment.  The  last  flight 
load  in  each  segment  is  assumed  to  be  a  minimum. 
As  an  option,  the  time  spent  in  a  given  flight 
may  be  selected  randomly.  A  selection  may  be 
made  for  each  mission  Interval  (or  weight 
interval)  in  a  mission. 

(d)  The  landing  Impact  follows  the  last  minimum 
flight  load.  This  is  a  one  time  loading  that 
Is  derived  from  a  distribution  of  sink  speeds, 
attitude,  etc. 

(e)  The  final  part  of  the  loading  sequence  is  the 
landing  taxi.  This  is  a  cyclic  loading  with 
the  first  load  a  maximum  followed  by  a  minimum 
and  then  another  maximum  etc.  It  Is  assumed 
that  the  last  load  Is  a  minimum.  The  time 
spent  In  landing  taxi  Is  that  which  is  appro¬ 
priate  to  the  mission  specified. 
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For  each  flight  load  occurrence  In  a  mission  or  weight  Interval 
there  will  be  an  associated  temperature  derived  from  the  velocity 
and  altitude  variations  within  that  segment.  It  will  be  assumed 
that  the  population  of  temperatures  is  Independent  of  the  population 
of  loads.  Consequently,  the  sampling  of  che  temperature  population 
does  not  depend  on  the  loading  that  ha.,  been  selected  for  that  time 
in  the  mission. 

The  following  definitions  describe  the  functions  required  to 
Implement  this  process. 

Suppose  each  of  a,  b,  Np,  and  Nm  a  positive  integer  such  that  a 
is  In  [1,  Np]  and  b  Is  In  [1,  Nm].  Further,  suppose  that  oab  is  the 

stress  at  the  ath  control  point  for  the  aircraft  resting  on**  its 
landing  gear  ready  to  fly  the  bth  mission. 

Suppose  that  My^  Is  a  positive  Integer  and  is  a  simple 

P  ,  p 

graph  such  that  the  x-projection  of  o^y  is  the  set  of  Integers  in 

[1.  ]  and  if  i  and  1*1  are  In 

P 


[1 »  Myy  ]  then  <4^  m  is 


TT  (1) 

K  r 

less  than  o*y  (j+|)*  is  a  finite  number  sequence  with  x- 


[1,  M#p  ]  and  nyyp  is  a  simple  graph  such  that  the  point 


ab 

TT 

.abl 


projection  [1,  M^yp].  How  suppose  that  1  is  a  positive  integer  in 


ab 


nT?p  (°TTp<i ) ) )  belon9S  to  nj?p  only  if  (o^p(D)  Is  the  number 

of  take-off  taxi  load  maximums  during  one  lifetime  at  the  stress 
oyyp(1)  for  the  ath  control  point  of  the  aircraft  flying  the  bth 
mission.  The  total  number  of  these  take-off  taxi  load  maximums  is 


.abl 


abl  /  ab 


»  rTTp  „abl  /  ab 

TT.  1  nTTp  <°TTp(1)) 


.abl 


PT 


1  -  1 


Therefore  there  are  n^y  selection  candidates  to  obtain  the 

PT 

the  first  one  of  the  M^yp  stresses.  A  selection  is  made  from  the 

population  of  Integers  In  [1,  n*y  ]  on  the  basis  that  If  1  and  j 

PT 


are  Integers  in  [1,  ]  then  the  probability  of  choosing  i  is 

PT 

equal  to  the  probability  of  choosing  j.  Now,  suppose  that  n^p  is 

a  simple  graph  with  x-projection  zero  plus  the  set  of  integers  in 
[1.  Myj  ]•  Also,  nf^J(o)  -  0  and  if  i  is  in  [1,  M^p]  then 


"SJm  -  l 


1 


nTTp  ^aTTp^k^* 


Therefore,  if  j  is  an  Integer  chosen  from  the  Integers  in  [1,  n2y^  ] 

Pf 


abl , 


then  there  exists  a  number  m  such  that  n^C«-l> 


,  t  <  abl 
3  *  nTTp(m). 


It  is  desired 


abl 


to  sample  the  population  of  Integers  in  [1,  n^fp  ] 


without  replacement.  Therefore,  before  the  second  selection  of  a  maximum 
stress  is  made  the  number  of  maximums  at  the  stress  o*yp(m)  must  be 
reduced  by  one.  Based  on  this  reduction,  the  simple  graphs  n^p 

and  n^p  are  transformed  into  n^p  and  n^p.  It  follows  then  if 


nTT  "  nTT  "1  an<l  J  Is  *n  integer  chosen  from  the  integers 
PT  PT 

in  [1,  n^p]  then  there  exists  a  number  m  such  that  njjp(m-l)  <  j  s  n^jp(m). 

This  process  can  be  repeated  until  all  of  the  maximums  have  been  selected. 
Suppose  that  Mj^n  Is  a  positive  Integer  and  o^n  Is  a  finite  number 

sequence  with  x-projectlon  [1,  M^].  Now  suppose  that  1  Is  a  positive 

Integer  in  [1,  M^n]  and  n^  is  a  simple  graph  such  that  the  point 
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(°TTnW’  nffj  (affn(1)))  belongs  to  nffj  only  If  nffj  (aj}n(1))  Is  the 


number  of  take-off  taxi  load  mlnlmums  during  one  lifetime  at  the  stress 

off  (1)  for  the  ath  control  point  of  the  aircraft  flying  the  bth  mission. 

The  total  number  of  these  take-off  taxi  load  mlnlmums  is 
Mab 

_abl  _  ™TTn  abl  /  ab  tt\\  .  „abl 


n 


n. 


n. 


=  l 
i 


_aDi  tao  ti\\  _  aD i 

,  nn  '“W1"  "  "tt 

I  n  PT 


In  this  case  there  are  the  same  number  of  selection  candidates  as 
for  the  take-off  taxi  load  maxlmums.  The  selection  from  the  mlnlmums 
Is  made  on  the  same  basis  as  from  the  maxlmums.  Suppose  that 

is  a  simple  graph  with  x-projectlon  zero  plus  the  set  of  Integers  In 
[1,  Further,  suppose  that  nff^(o)  *  0  and  if  1  Is  in  [1» 


then 


abl / , , 
nTTn(i) 


1 


nTTn  (°TTn(k^ 


Therefore,  If  j  is  an  Integer  chosen  from  the  Integers  In 
[1.  nff1  ]  then  there  exists  a  number  m  such  that  <  j  -  ^jT^m) 

fly  ■»  *1 


As  In  the  case  of  the  maxlmums  It  is  desired  to  sample  the  mlnlmums 

without  replacement.  Consequently,  the  simple  graphs  nff^  and  are 

n  n 

transformed  to  nif2  and  nff2.  etc. 

n  n 

Suppose  that  each  of  c  and  Is  a  positive  Integer  such  that  c  is 

In  [1,  h£].  Further,  suppose  that  M?bc  Is  a  positive  Integer  and 

P 

o*bc  Is  a  finite  number  sequence  with  x-projectlon  [1,  Mpbc].  flow. 

P  P 

suppose  that  1  Is  a  positive  integer  In  [1,  HSbc]  and  n?bc^  Is  a 

P  P 


6 


>i  in  '  . 


simple  graph  such  that  the  point  (apb(.j,  npbcl(<Jpb(j)))  belongs  to 

npbc^  only  If  njibc^{opb9. »)  Is  the  number  of  flight  "positive  g  loads" 
rP  P  p'1' 

during  one  lifetime  at  the  stress  op  for  the  ath  control  point  of 

the  airplane  flying  the  cth  interval  of  the  bth  mission.  The  term 

"positive  g  loads"  refers  to  those  maximum  load  occurrences  corresponding 

to  nz  In  excess  of  one. 

The  total  number  of  these  positive  g  loads  Is 
Mabc 

rf*1  .  IFp  n»bcl  („"b<:  ) 

FPT  1*1  *P  Fp(1) 

For  the  first  load  selection  after  take-off  there  are  n^1 


selection  candidates  to  obtain  one  of  Mp  stresses.  A  selection 

P 

is  made  from  the  population  of  Integers  In  [1,  n2b11]  on  an  equal 

PpT 

probability  basis  (as  for  the  taxi).  Now,  suppose  that  npbc^  Is  a 

P 

simple  graph  with  x-projection  zero  plus  the  set  of  integers  In 
[1,  Mjbc].  Also,  npb/J)  -  o  and  if  1  is  In  [1,  Mj?*]  then 


«!bc'  .  I  n*bc’  (.f^k)) 

Fp(i)  k  ■  1  Fp  Fp 


P(i)  k  *  l  p  p 

Thus,  If  j  Is  an  Integer  chosen  from  the  Integers  In  [1,  n?b^] 

PT 

then  there  exists  a  number  m  such  that  <  j  i 


7 


This  population  Is  also  sampled  segment  by  segment  without 
replacement  until  all  of  the  maximums  are  selected. 


Suppose  that  M!bc  is  a  positive  Integer  and  opbc  Is  a  finite 
n  rn 


number  sequence  with  x-projectlon  [1,  M?bc  ].  Further,  suppose  that 

n 

1  Is  a  positive  integer  In  [1,  Mfbc]  and  njbc1  is  a  simple  graph  such 

n  n 

that  the  point  (a?b/n,  njbcl  (o?b?n))  belongs  to  n2bc1  only  if 
n'  ;  n  n'i;  rn 

n?bc1 (cpb9 . . )  is  the  number  of  "negative  g  loads"  during  one  lifetime 
n  n*  ' 


at  the  stress  apb^j  for  the  ath  control  point  of  the  airplane  flying 


the  cth  interval  of  the  bth  mission.  The  term  "negative  g  loads"  refers 
to  those  minimum  load  occurrences  corresponding  to  nz  less  than  one. 


The  total  number  of  these  negative  g  loads  is 


M! 


,abc 


^abd  .  j  n  abcl  (abc  (1)) 

hnT  1-1  *'n  Fn 


It  Is  assumed  that  each  positive  g  loading  is  followed  by  either 
a  negative  g  load  or  by  a  "one  g  load"  (a  "one  g  load"  Is  a  load  cor¬ 
responding  to  an  nz  minimum  equal  to  one). 


The  total  number  of  one  g  loads  for  the  ath  control  point  of  the 

aircraft  flying  the  cth  Interval  of  the  bth  mission  is 

_abcl  .  „abcl  „abcl 
"F.  nF„  -  nF. 


S 


The  stress  spectrum  for  the  one  g  loads  will  be  derived  from  the 
velocity,  altitude,  and  weight  interval  load  occurrences  that  were 
used  for  the  positive  g  load  stress  spectrum.  It  is  assumed  that 
for  each  load  occurrence  the  normal  load  factor  will  be  equal  to  one. 
To  obtain  the  correct  nunber  of  minimums  from  this  spectrum,  the 
number  of  loadings  at  each  stress  level  will  have  to  be  multiplied 
by  the  ratio  nabc1/np  c  .  The  resultant  one  g  maneuver  load 

FgT  Pl 

spectrum  can  be  defined  as  follows: 

Suppose  MpgC  is  a  positive  integer  and  o*gC  is  a  finite 
number  sequence  with  x-projectlon  the  interval  [1,  Also, 

suppose  that  i  is  an  integer  in  [1,  M^c]  and  npgCl  is  a  simple 
graph  such  that  the  point  (opgC(1),  npgCl(<JpgC  (i ) ))  belongs  to 

nFgCl  only  nFgCl^°FgC^^  the  number  "one  9  loads" 

during  one  lifetime  at  the  stress  crpgC(1)  for  the  ath  control 

point  of  the  airplane  flying  the  cth  interval  of  the  bth  mission. 

Note  that 


_abcl 

nF_ 


»2bc'  („*»(,» 
g  g 


The  selection  of  the  flight  load  minimum  is  complicated  by  the 
fact  that  it  may  come  from  one  of  the  negative  g  loads  or  it  may  come 
from  one  of  the  one  g  loads.  This  condition  may  be  handled  as  follows: 
Suppose  abc1  is  a  simple  graph  with  x-projection  zero  plus  the  set 

nFn 

of  integers  in  [1,  MpgC  +  M*^c  ].  If  i  is  in  [1,  MpgC]  then 


(1)  is  r  n*!*1  (o£c(k)).  If  i  is 


C^g0  +  1  •  MFgC+  ^  then  nfjc1(i)  is 


1  *  N 

+  i  njjc1  (oj5c(h)).  Also,  npJcl(o)  -  o.  Therefore, 


f 


for  the  first  flight  load  minimum  selection.  If  j  Is  an  Integer 
chosen  from  the  Integers  in  [1,  abl]  ]  then  there  exists  a  number 

nppT 

m  such  that  hpn^(m-l)  <  j  <  nFb^m).  Again,  the  sampling  is  made 
without  replacement  so  that  the  number  of  mlnlmions  is  reduced  by 


one.  The  transformations  of  n*^1,  nfb11,  and 


+a  nab12  nab]2 

to  nFg  ,  nFn  , 


and  follows  from  this  reqirement.  This  population  Is  sampled 

segment  by  segment  until  all  of  the  mini mums  are  drawn. 

The  landing  Impact  stress  spectrum  will  be  sampled  once  per 
flight  without  replacement.  This  spectrum  can  be  defined  as  follows: 

Suppose  Mj^  Is  a  positive  number  and  ofb  is  a  finite  number  sequence 

with  x-projectlon  [1,  Mfb  ].  Further,  suppose  j  is  a  positive  integer 
In  [1,  Mfb  ]  and  nfb^  is  a  simple  graph  such  that  the  point  (ofb  (j), 

"lJ1  (olJ  belongs  to  n*bl  only  If  n*bl  (ofj  (j))  is  the  number 
of  landing  impact  minlmums  during  one  lifetime  at  the  stress 
o^(j)  for  the  ath  control  point  of  the  airplane  flying  the  bth 
mission. 

Suppose  that  Mfb  is  a  positive  Integer  and  o?b  is  a  finite 

P  "  p 

number  sequence  with  x-projectlon  [1,  ].  Further,  suppose  that 


1  is  a  positive  Integer  in  [1,  M*b  ]  and  nft1  Is  a  simple  graph  such 

P  p 

that  the  point  (of®  (1),  n^1  (oLT  (1)))  belongs  to  n*b  only  if 

p  p  p  ‘■'p 

nLTp  (°LT  *be  nunt>er  of  landing  taxi  load  maxlmums  during 

one  lifetime  at  the  stress  of!j  (i)  for  the  ath  control  point  of  the 

P 
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aircraft  flying  the  bth  mission.  The  total  number  of  these 
maxi mums  Is 

yab 

n,  T 

„abl  ,  _  p  nabl  ,  ab  ,*u 

nLT  *  "  (oLTb{<))- 

PT  11  LTp  P 

The  landing  taxi  stress  selection  Is  analogous  to  that 
selection  used  for  the  take-off  stresses.  In  this  case  there 

are  nfE1  selection  candidates  to  obtain  the  first  one  of  the 


stresses.  This  selection  Is  made  on  an  equal  probability 
P  aK1 

basis.  Now  suppose  that  n^y  Is  a  simple  graph  with  x-projectlon 

zero  plus  the  set  of  Integers  In  [1,  M?y  ].  Also,  nfy^(o)  *  0  and 

k  L  P  L,p 

If  1  Is  in  [1,  M?S  ]  then 


nff’O)  -  l  nft1  (<$  (k)) 

LTp  k-1  LTp  LTp 


ah  I 

The  simple  graph  nLT  Is  used  to  determine  the  stress  that 
P  .m 

corresponds  to  the  selection  from  the  njj  Integers. 

PT 

Again,  the  sampling  is  done  without  replacement  until  the 
entire  population  has  been  selected. 

Suppose  M*y  Is  a  positive  integer  and  o*y  Is  a  finite 
n  .u  n 

number  sequence  with  x-projectlon  [1,  M|j  ].  Further,  suppose 

that  1  Is  a  positive  Integer  In  [1,  m£j  ]  and  n^y1  Is  a  simple 

n  n 

graph  such  that  the  point  (ojj  (1),  n*y*  (o*y  (1)))  belongs  to 

n  n  n 

njy1  only  If  n*y  (o*y  (1))  Is  the  number  of  landing  taxi  load 
n  n  n 


<4f4’.v 


mlnlmuras  during  one  lifetime  at  the  stress  (1)  for  the  ath 

L  n 

control  point  of  the  aircraft  flying  the  bth  mission.  The  total 
number  of  these  maxi mums  Is 

»ab 

abl  «LT  „abl  <  /,\\ 

LTn  lt  {°lt  (i)) 

nj  i=l  n  n 

The  minimum  stress  selection  Is  made  on  the  same  basis  as 

the  maximum  stress  selection.  For  this  purpose  suppose  that 

is  a  simple  graph  with  x-projection  zero  plus  the  set  of  integers 

In  [1,  M*b  ].  Also,  nj|bl(o)  =  0  and  if  1  is  in  [1,  M®b  ]  then 
n  n  n 


nLT1(i>  *  *  "ff1  (°LT  <k» 
L  n  k=1  L  n  L  n 


The  sampling  for  the  minimum  stresses  is  done  without  re¬ 
placement  until  all  of  the  candidates  have  been  selected. 

To  include  the  temperature  effects  on  the  test,  the  following 
functions  need  to  be  defined.  Suppose  that  ^abc  Is  a  positive 

integer  and  T^bc  is  a  finite  number  sequence  with  x-projectlon 

[1,  M*bc],  Further,  suppose  that  1  Is  a  positive  Integer  In 

[1,  M*bc]  and  n*bcl  Is  a  simple  graph  such  that  the  point 

(lpbc(1),  n*bcl(Tjbc(1)))  belongs  to  nfbc1  only  If  n*bcl  (T*bc(i)) 
Is  the  number  of  occurrences  during  one  lifetime  of  the  temperature 
lpbc(1 )  at  the  ath  control  point  of  the  airplane  flying  the  cth 
Interval  of  the  bth  mission. 

The  total  number  of  these  temperature  occurrences  must  agree 
with  the  number  of  load  occurrences.  That  Is, 
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Isabel  .  tF  „abd  (Tabc(1))  .  2„abcl 
1  i-i  hr 

The  temperature  sampling  is  made  without  replacement  first 
from  the  set  of  integers  in  [1,  n|b^].  For  this  purpose  suppose 

r7 

npbc1  is  a  simple  graph  with  x-projection  zero  plus  the  set  of 
integers  in  [1,  Mj!bcJ.  Also,  npbc^(o)  =  0  and  if  i  is  in  [1,  Mabc] 


npbc1(i)  =  Z  nabcl  (Tpbc(k)) 
h  k=l  F 

Therefore,  if  j  is  an  integer  chosen  from  the  integers  in 

H  nab11 

L  ’  F  ]  then  there  exists  a  number  m  such  that 
T 

nab11  (m-1 )  <j<njbll(m) 


This  population  is  sampled  segment  by  segment  until  all 
temperatures  have  been  selected. 

The  special  case  of  the  touch-and-go  landing  must  be  considered. 
For  this  case  it  is  assumed  that  there  are  two  loads.  First,  there 
is  a  landing  impact  load  followed  by  a  "one  g"  load.  The  first  touch- 
and-go  in  a  given  mission  is  assumed  to  be  preceded  by  a  flight  load 
minimum.  If  there  are  repeated  touch-and-go  landings  it  is  assumed 
that  there  are  no  flight  loads  between  them. 

The  required  functions  for  this  case  are  defined  as  follows: 
Suppose  uab  is  a  positive  Integer  and  ab  is  a  finite  number 
mtg1  °TG. 

sequence  with  x-projection  [1,  M^|?  ].  Further  suppoie  that  j  Is  a 

positive  Integer  in  [1,  ]  and  naj^  is  a  simple  graph  such  that 

the  point  (aaj^(j),  njg|  (o^U)))  belongs  to  njj?j  only  If  n^J 

(°TG  the  number  of  touch-and-go  landing  Impact  loads  during 

*  ah  . 

one  lifetime  at  the  stress  oyj^(j)  for  the  ath  control  point  of  the 
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aircraft  flying  the  bth  mission.  The  total  number  of  these  landing 
impact  loads  Is 

M?b 

abl  ^Gi  abl  ,  ab  , . u 

nTG,  "  1  nTG.  ^°Tg/j^ 

’T  j  =1  1  1 

Now,  suppose  that  hjq1  is  a  simple  graph  such  that  the  x-projectlon 
of  i>j^  is  zero  plus  the  integers  in  [1,  ].  Also,  nj^(o)  *  0  and 

if  j  is  in  [1 ,  Mjg  ]  then 


=  L  nTG^ 


i 


k*l 


The  sampling  of  the  population  of  integers  in  [1,  nyj?1] 

T 

without  replacement  as  In  the  case  of  the  previously  defined  normal 
landing  impact. 

For  the  following  one  g  load  suppose  that  MyjJ  is  a  positive 

ab  9 

integer  and  a?-  Is  a  finite  number  sequence  with  x-projection 

9 

[1,  M^q  ].  Further,  suppose  that  1  is  a  positive  integer  In  [1,  ] 

9  9 

and  n?-1  Is  a  simple  graph  such  that  the  point  (oyG  (i),  nJg1  (oTG  (1))) 
g  9  9  9 

belongs  to  nyjJ1  only  If  nyj?1  (a™  (1))  Is  the  number  of  one  g  load  con- 
9  9  9 

dltions  (after  a  touch-and-go  landing)  during  one  lifetime  at  the  stress 

°TG  (0  for  *th  control  point  of  the  aircraft  flying  the  bth  mission. 
9 

The  total  number  of  these  loads  must  be  the  same  as  the  touch-and-go 
impact  loads.  Therefore 


.abl 

TG_ 


■  t  *  »•£'<„•]!  (0)  <  nIb’ 


i-1 


TG, 
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To  fascillitate  the  sampling  process  suppose  that  is  a 

9 

simple  graph  such  that  the  x-projectlon  of  is  zero  plus  the 

9 

set  of  integers  in  [1,  M^!  ].  Further  njj^(o)  ■  0  and  if  i  is  in 

9  9 

tl.  ]  then 
9 

nTG^^  *  l  ,  nTG^  ^°TG 
g  k=l  g  g 

The  sampling  of  the  population  of  integers  in  [1,  n—1  ]  is  made 

9 

without  replacement. 


SECTION  III 


DESCRIPTION  OF  COMPUTER  PROGRAM 


I  NOTATION 

NCASE  =1  if  flight  maneuver  stress  spectra  Is  derived  from 
recorded  data 

=  2  if  flight  maneuver  stress  spectra  is  derived  from 
mission  profiles 

NM  -  Number  of  missions 

NRANM  =1  if  missions  are  randomly  selected 

=2  if  mission  selection  Is  nonrandom 

NMTMS  =1  if  the  number  of  maximum  flight  load  occurrences  per 
mission  is  nonrandom.  If  NMTMS  >1  then  NMTMS 
is  the  number  of  possible  candidates  for  random 
selection  of  the  number  of  maximum  flight  load  occurrences 
per  mission 

NFLT  -  Number  of  flights  in  one  lifetime 

NHRS  -  Number  of  flight  hours  in  one  lifetime 

MISRPT  -  Number  of  times  a  mission  sequence  occurs  in  a  lifetime 
(for  nonrandom  mission  selection) 

MNUMS  -  Number  of  different  mission  selections  in  a  sequence 
(for  nonrandom  mission  selection) 

NS ( 1 )  -  Number  of  mission  Intervals  in  the  ith  mission 
(excluding  take-off  and  landing) 

NTG(1)  -  Number  of  touch-and-go  landings  in  the  1th  mission 

NTEMP(i)  *  1  if  a  temperature  profile  is  not  included  In  the  1th 
mission 

=2  if  a  temperature  profile  is  Included  In  the  1th  mission 

NTTF(1)  -  Number  of  maximum  take-off  taxi  load  occurrences  per  flight 
for  the  1th  mission 

NFF(k,1,j)  -  Number  of  maximum  flight  load  occurrences  per  flight  for 
the  kth  selection  from  NMTMS  possible  candidates  for 
the  jth  Interval  of  the  1th  mission 

NITF(I)  -  Number  of  maximum  landing  taxi  load  occurrences  per  flight 
for  the  1th  mission 
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MNUM(i)  -  The  ith  mission  number  selected  (non randomly)  in  amission 
sequence  that  occurs  MISRPT  times  in  a  lifetime 

NUMRPT(i)  -  The  number  of  entries  of  MNUM(i)  within  a  mission  sequence 

NFRAC(k,  i,  j)  -  Number  of  times  per  lifetime  that  the  kth  selection  from 
NMTMS  possible  candidates  for  the  jth  Interval  of  the  1th 
mission  that  the  number  NFF(k,  1,  j)  is  selected 

NUMML(i)  -  Number  of  times  the  1th  mission  occurs  in  one  lifetime 

MTTP(i)  -  Number  of  maximum  take-off  taxi  load  Intervals  for  the  ith 
mission 

MTTN(i)  -  Number  of  minimum  take-off  taxi  load  Intervals  for  the  1th 
mission 

MFP(1,  j)  -  Number  of  positive  g  load  intervals  for  the  jth  interval  of 

the  ith  mission 

MFN(1 ,  j)  -  Number  of  negative  g  load  Intervals  for  the  jth  interval 
of  the  1th  mission 

MFG(1,  j)  -  Number  of  one  g  load  Intervals  for  the  jth  Interval  of 
the  1th  mission 

MLI(1)  -  Number  of  landing  Impact  load  Intervals  for  the  1th  mission 

MLTP ( 1 )  -  Number  of  maximum  landing  taxi  load  intervals  for  the  1th 
mission 

MLTN(1)  -  Number  of  minimum  landing  taxi  load  Intervals  for  the  1th 
mission 

MF(1 ,  j)  -  Number  of  temperature  intervals  for  the  jth  Interval  of 
the  1th  mission 

MTGI(1)  -  Number  of  touch-and-go  landing  Impact  load  Intervals  for 
the  1th  mission 

MTGG(1)  -  Number  of  touch-and-go  landing  one  g  load  Intervals  for  the 
1th  mission 

HTTP(1,  j)  -  Number  of  maximum  take-off  taxi  load  occurrences  In  the 
jth  load  Interval  for  the  ith  mission 

NTTN(1,  j)  -  Number  of  minimum  take-off  taxi  load  occurrences  In  the 
jth  load  Interval  for  the  1th  mission 

NFP(1,  j,  k)  -  Number  of  positive  g  load  occurrences  In  the  kth  load 
Interval  for  the  jth  Interval  of  the  ith  mission 


NFN(1,  j,  k)  -  Number  of  negative  g  load  occurrences  In  the  kth 

load  Interval  for  the  jth  Interval  of  the  1th  mission 

NFG(1 ,  j,  k)  -  Number  of  one  g  load  occurrences  In  the  kth  load  Interval 
for  the  jth  Interval  of  the  1th  mission 

NLI{1,  j)  -  Number  of  landing  Impact  load  occurrences  In  the  jth  load 
Interval  for  the  1th  mission 

NLTP{1,  j)  -  Number  of  maximum  landing  taxi  load  occurrences  in  the 
jth  load  Interval  for  the  1th  mission 

NLTN(1,  j)  -  Number  of  minimum  landing  taxi  load  occurrences  in  the 
jth  load  Interval  for  the  1th  mission 

NF (1 ,  j,  k)  -  Number  of  temperature  occurrences  In  the  kth  temperature 
Interval  for  the  jth  Interval  of  the  1th  mission 

NTGI(1 ,  j)  -  Number  of  touch-and-go  landing  Impact  load  occurrences  per 
lifetime  In  the  jth  load  Interval  for  the  1th  mission 

NTGG(1 ,  j)  -  Number  of  touch-and-go  landing  one  g  load  occurrences  per 
lifetime  In  the  jth  load  Interval  for  the  1th  mission 

TL  -  Temperature  for  all  ground  load  conditions 

SR(i)  -  Stress  for  the  take-off  configuration  for  the  1th  mission 
(aircraft  resting  on  the  landing  gear) 

STTP(1 ,  j)  -  Stress  in  the  jth  maximum  take-off  taxi  load  Interval 
for  the  1th  mission 

STTN(l.j)  -  Stress  In  the  jth  minimum  take-off  taxi  load  Interval 
for  the  1th  mission 

SFP(1,  j,  k)  -  Stress  In  the  kth  positive  g  flight  load  Interval  for 
the  jth  Interval  of  the  1th  mission 

SFN(1,  j,  k)  -  Stress  In  the  kth  negative  g  flight  load  Interval  for 
the  jth  Interval  of  the  1th  mission 

SFG(1,  j,  k)  -  Stress  In  the  kth  one  g  flight  load  Interval  for  the 
jth  interval  of  the  1th  mission 

SLI (1  •  j)  -  Stress  In  the  jth  landing  impact  load  Interval  for  the 
ith  mission 

SLTP(1,  j)  -  Stress  In  the  jth  maximum  landing  taxi  load  Interval 
for  the  1th  mission 

SLTN(1,  j)  -  Stress  In  the  jth  minimum  landing  taxi  load  Interval 
for  the  1th  mission 


TF{1,  j,  k)  -  Temperature  In  the  kth  flight  temperature  Interval 
for  the  jth  Interval  of  the  1th  mission 

STGI(1,  j)  -  Stress  In  the  jth  touch-and-go  landing  Impact  load 
Interval  for  the  1th  mission 

STGG(1,  j)  -  Stress  In  the  jth  touch-and-go  one  g  load  Interval 
for  the  1th  mission 


2  COMPUTER  FLOW  DIAGRAM  AND  PROGRAMS 

The  computer  routine  was  coded  In  FORTRAN  Extended  Language 
with  the  main  program  and  subroutines  arranged  as  follows: 

SPECF  Program  (See  Appendix  A) 


MAIN  -  Main  Program  -  Sets  NZERO  to  zero  and  transfers  control  to  GUIDE 

GUIDE  -  Subroutine  -  Initially  zeros  Input  and  output  numbers  and  after 
first  case  zeros  output  numbers  before  the  calculations  are 
performed  GUIDE  Is  the  subroutine  for  transferring  control 
to  INPUT,  COW,  and  CALC  In  turn. 

INPUT  -  Subroutine  -  Reads  In  all  Input  data.  The  details  of  the  data 
Input  are  discussed  later  In  this  section 

INAID  -  Subroutine  -  Called  by  INPUT  to  write  certain  input  data. 

INAID  sets  NZERO  -  1  for  control  In  GUIDE 


PAGEHD  -  Subroutine  -  Writes  out  page  heading  Including  run  Identi¬ 
fication,  date,  and  page  number 


COMP  -  Subroutine  -  Checks  the  Input  data  for  compatibility  as  follows: 


For  1  in  [1 ,  NM] 
NTTF(1)  •  NUMML(I) 


MTTP(1) 

I  NTTP(1,  j) 

J-l 
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MTTN(1) 

NTTF(1 )  •  NUWL(i)  -  t  NTTN(i,  j) 

ML 1(1) 

NUMML(I)  «  £  NLI ( 1 ,  j) 

j»1 

MLTP(1) 

NLTF(1)  •  NUMML ( 1 )  «  £  NLTP(1,  j) 

j»1 

ML  1(1) 

NLTF(1)  *  Nim(i)  *  £ 

MTGI(1) 

NTG( 1 )  •  NUMML(i)  *  £  NTGI(1,  j) 
j-1 

MTGG(1) 

NTG(1)  •  NUMML ( 1 )  -  £  NTGG(1,  j) 

j“1 

NM 

NFLT  -  £  NUMML  ( 1 ) 

1-1 


For  j  In  [1,  NS (1)3 

NMTMS  MFP{1 ,  j) 

£  NFF(1,  j,  k)  *  NFRAC  (i,  j,  k)  »  £  NFP(1,  j.  k) 

k«l  k-1 

NMTMS  MFN{1 ,  j)  MFG(l.j) 

£  NFF(1,  j,  k)  •  NFRAC(1,  j,  k)  -  £  NFN(1,j,  k)  +  £  NFG(1,j,k) 

k-1  k-1  k»l 


If  any  of  the  above  conditions  are  violated, an  error  signal  will  be 

printed  and  GUIDE  will  be  called. 

CALC  -  Subroutine  -  Computes  and  prints  fllght-by- flight  load  sequence 

SEARC1  -  Subroutine  -  Makes  random  selection  of  stress  where  the 
number  of  occurrences  per  lifetime  Is  a  one  dimensional 
array 

SEARC2  -  Subroutine  -  Makes  random  selection  of  stress  where  the 
number  of  occurrences  per  lifetime  Is  a  two  dimensional 
array 

SEARC3  -  Subroutine  -  Makes  random  selection  of  stress  of  temperature 
where  the  number  of  occurrences  per  lifetime  Is  a  three 
dimensional  array 
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3  EQUIVALENCE  TABLES 

All  Input  numbers  for  this  routine  are  placed  In  blank  common. 
All  Input  floating  point  numbers  are  called  parameters  and  are  con¬ 
tained  In  P  (dimensioned  13.000).  All  Input  fixed  point  numbers  are 
called  Integers  and  are  contained  In  N  (dimensioned  13.500).  To 
make  the  program  usable,  EQUIVALENCE  statements  are  used  to  give  the 
P  and  N  numbers  more  recognizable  names.  The  SPECF  program  parameter 
and  Integer  tables  are  given  below. 


EQUIVALENCE  TABLE  -  P  (Parameters) 


p 

Dim 

Term 

P 

Dim 

Term 

1 

(1) 

TL 

101 

111 

115 

354 

415 

654 

715 

3114 

3115 

5514 

5515 

7914 

7915 
8154 
8215 
8454 
8515 
8754 
8815 

11214 

11215 
11454 
11515 
11754 
11815 

12114 

12115 
12414 

(12) 

(12.20) 

(12,20) 

(12,  10,  20) 

(12,  10,  20) 

{12,  10,  20) 

(12.  20) 

(12.  20) 

(12,  20) 

(12.  10,  20) 

(12.  20) 

(12.  20) 

(300) 

(300) 

SR(1) 

SR(12) 

STTP  (1,1) 
STTP(12,  20) 
STTN  (1,1) 
STTN(12,  20) 
SFPtl,  1.  lJ 
SFP  12,  10,20) 
SFN1,  1,  1) 
SFN  12,  10,20) 
SFGtl,  1.  1} 
SFG(12,  10.20) 
SLI(1,  1) 

SLI  12.  20) 
SLTP  1,  1) 
SLTpne,  20) 
SLTNtl ,  ,1) 
SLTN(12,  20) 

TF  1.  1.  1) 
TF(12,  10,20) 
STGI(1,  1) 

STGI  12,  20) 

STGGil,  1) 

STGG(12,  20) 

STRESS(l) 

STRESS(300) 

TEMP(l) 

TEHP(300) 

EQUIVALENCE  TABLE  -  N  (Integer) 


I  DENT 

NPF1 

NPF2 

NPF3 

NPF4 

NPF5 

NPF6 

MONTH 

DAY 

YEAR 

NCASE 

NM 

NRANM 

NMORE 

NMTMS 

NFLT 

NHRS 

MISRPT 

MNUMS 


NZERO 

NPAGE 


N 

Dim 

115 

354 

(12.  20) 

415 

654 

(12,  20) 

715 

3114 

(12,  10.  20) 

3115 

5514 

(12,  10,  20) 

5515 

7914 

(12.  10,  20) 

7915 

8154 

(12.  20) 

8215 

8454 

(12,  20) 

8515 

8754 

(12,  20) 

8815 

11214 

(12,  10.  20) 

11215 

11454 

(12,  20) 

11515 

11754 

(12.  20) 

11815 

11826 

(12) 

11830 

11841 

(12) 

11845 

11964 

(12,  10) 

11965 

12084 

(12,  10) 

12085 

12204 

(12.  10) 

12205 

12216 

02) 

12220 

12231 

02) 

NTTP(1,  1) 
NTTP(12,  20) 
NTTNlI ,  1) 
NTTN(1,2,  20) 
NFP  (1,  1,  1) 
NFP(12,10,  20) 
NFNfl ,  1.  1) 
NFN(12,  10,20) 
NFG(1 ,  1.  1) 
NFG(12,?Q.  20) 
NLljl,  1) 
NLI(12,  20) 
NLTP(1,  1) 
NLTP(12,  20) 
NLTN(1,  1) 
NLTN(12,  20) 
NFM.  1,  1) 
NF(12»  10.20) 
NTGI  1,  1) 
NTG1  12.  20) 
NTGG  1.  1) 
NTGGi 12,  20) 
MTTPil ) 
MTTP(12) 
MTTN(l) 
MTTN(12) 
MFP(1,  1) 
MFP(12,  10) 
MFN  1,  1) 
MFN(12,  10) 
MFG  1,  1) 

MFG( 12 ,  10) 
MLI(1) 

MLI(12) 

MLTP(l) 

MLTP(12) 


-■'saam! ■  -  --- 


EQUIVALENCE  TABLE  -  N  - 


N 

Din 

Term 

12235 

(12) 

MLTN(I) 

12246 

HLTN(12) 

12250 

(12.  20) 

MF  1.  1)  , 

12369 

HF(12.  10) 

12370 

(12) 

MTGI(l) 

12381 

MTGl(l2) 

12385 

(12) 

NTGGcl ) 

12396 

MTGG(12) 

12400 

12411 

(12) 

N S(1) 

NS(12) 

12415 

(12) 

NTG(l) 

12426 

NTG(12) 

12430 

12441 

(12) 

NTEMP(I) 

NTEHP(12) 

12445 

12456 

(12) 

NTTF(1 ) 

NTFF(12) 

12460 

12819 

(*.  1?.  TO) 

NFF(1 ,  1,  1) 
NFF(3.  12.  10) 
NLTF(1) 

12820 

(12) 

12831 

NLTF(12) 

12835 

(50) 

HNUM(l) 

12884 

MNUM(50) 

12885 

(50) 

NUHRPT(1 ) 

12934 

NUMRPT(50) 

12935 

(3.  12.  10) 

) 

NFRACM,  1.1) 

13294 

13295 

NFRACC3,  12,  10) 
NUMML(l) 

(J2)  / 

13306 

NUWL(12) 

(Intege 


4  INPUT  DATA 

All  of  the  Input  data  described  below  Is  read  Into  the  program 
by  means  of  the  subroutine  INPUT.  This  program  Is  a  general  purpose 
routine  to  read  the  P  (parameters)  and  N  Integers.  There  are  several 
options  by  which  this  may  be  done  by  this  routine.  A  suggested  deck 
arrangement  Is  give*  as  follows: 

1415  Format 


I  DENT 

1 

1 

NPF3 

1  0 

NPF5 

NPF6 

MONTH 

DAY 

YEAR 

NCASE 

NM 

NRANM 

0 

515  Format 


NMTMS 

NFLT 

NHRS 

MISRPT 

MNUMS 

72H  Format 
Run  Description 


72H  Format 
Run  Description 


315  Format 


101 


100  +  NM 


1 


6E10.3  Format 

Q _ SR(1)  -  SR(NM) _ 


15  El 5. 7  Format 


315  Format 


1 21 5  Format 


NS(1 )  -  NS(NM) _ 

( NS ( 1 )  must  not  exceed  10) 


315  Format 


12415 


12414  +  NM  1 


1215  Format 
NT6(1)  -  NTG(NM)  ~~ 


315  Format 


|  12430 


12429  +  NM 


1215  Format 

1  NTEMP(l)  -  NTEMP(NM) 


315  Format 


12445 


12444  +  NM 


1215  Format 
NTTF(l)  -  NTTF(NM) _ 


315  Format 

12820  |  12819  +  NM  1  7 


1215  Format 
NLTF(l)  -  NLTF(NM) 


315  Format 


If  NRANM  -  1  go  to  (bb),  if  NRANM  -  2  go  to  (aa) 

315  Format 


(bb)  315  Format 


(MTTP(i)  must  not  exceed  20) 

315  Format 


I 


(MTTN( 1 )  must  not  exceed  20) 


315  Format 


12220  1  12219  +  NM  fT 


1215  Format 


MLTP(l)  -  MLTP(NM) 


315  Format 


112235  |  12234  +  NM  1  1  1 


1215  Format 


MLTN{1)  -  MLTN(NM) 


If  1  is  in  [1,  NM]  and  NTG(i)  *  0  go  to  (b),  otherwise  go  to  (a) 

(a)  315  Format 

j  12370  i  12369  +  NM  I  V] 

1215  Format 

MTGI(1 )  -  HTSI(NM)  1 

(MTGl(i)  must  not  exceed  20) 


27 


315  Format 


i 

i 

’» 


12385 


12384  +  NM 


1 


1215  Format 
MTGG(l)  -  MTGG(NM) 

(MTGG(1 )  must  not  exceed  20) 

15  Format 

(b)  \W\ 


715  Format 


115  1  1  MTTP(l)  1 

1  12 

i  20 

6E10.3  Format 

snp(i,  i)  -  snp(i,  mttp(D) 

• 

• 

• 

715  Format 

115  1  1  MTTP(i)  1 

12 

20  | 

6E10.3  Format 

STTP(t,  1)  -  STTP (1 ,  MTTP(l)) 


715  Format 


115 

□_ 

NH 

MTTP(NM)  j 

1 

12 

20  | 

6E10.3  Format 


STTP(NM,  1)  -  STTP(NM,  MTTP(NM)) 


15  Format 

m 


715  Format 


415 

1 

1 

MTTN(l) 

1 

12 

20 

6E10.3  Format 


STTN(1,  1)  -  STTN(1 ,  MTTN(1 )) 


715  Format 


415 

1 

i 

MTTN(  i  ) 

1 

12 

L_20j 

6E10.3  Format 


STTN(i,  1)  -  STTN(1,  MTTN(i)) 


715  Format 


415  1 


NM 


MTTN(NM)  |  1  \  12 


20 


6E10.3  Format 


STTN(NM,  1)  -  sm(NM,  MTTN(NM)) 


15  Format 


NM 

Z  NS(1 ) 
1«1 


715  Format 

715  |  1  [  1  MFP(1,  1)  |  12  |  10  |  20 


29 


6E10.3  Format 

SFP(1,  1,1)  -  SFP(1  >1 ,  MFP(1,  1)) 


715  Format 

715  1  i  1  j  1  MFP(  i  ,j)l  12  1  TO  1  20 


6E10.3  Format 

SFP(1,  j.  1)  -SFP(1,j,  MFP(t ,  j)) 


715  Format 

715  j  NM  1  NS(NM)  HFP(ffl,NS(NMa  12  1  10  t  20 


6E10.3  Format 

SFP(NM,  NS(NM),1)  -  SFP(NH,  NS(NM),  HFP(NM,  HS(NM)) 


15  Format 
"NM 

Z  NS(1 ) 
1-1 


715  Format 

3115  1  HU  ^(1,  1  12  |  10  1  20 

6E10.3  Format 

SFN (1 ,  1,1)  -  SFN(1 ,  1,  MFN~(1,  1)) 


715  Format 

3115  1  1  |  J  |  HFW(1,  J)  1  12  j  10  j  20 


30 


6E10.3  Format 


$FN(1,  j,  1)  -  SFN ( 1 ,  j,  MFN(1,  j)) 


715  Format 


I  3115 


NM 


NS(NM)  MFN(NM,  NS(NM))  |  12  |  10  1  20 1 


6E10.3  Format 

SFN(NM,  NS(NM) ,  1)  -  SFN(NM,  NS(NM),  MFN(NM,  NS(NM)) 


15  Format 


NM 

£  NS(i) 
1-1 


715  Format 


1  |  1  1  MFG(lt  1) 


5515 


12 


10 


20 


6E10.3  Format 

i"  SFG(1.  1.  1)  -  SFG(1,  1,  MFG(1 ,  1))~ 


> 


715  Format 


5515  1 


MFG(1,  j)  12  10  20 


§E10.3  Format 

SFG(i ,  j,  1)  -  $FG(1 ,  j,  MFG(i ,  j)) 


715  Format 


5515 


NM 


NS(NM) 


MFG(NM,  NS(NM))  12 


10  20 


i 

I 

i 

! 


6E10.3  Format 


SFG(NM,  MS(NM),  1)  -  SF6(»H,  HS(NH),  HF6(NH,  NS(NH)) 


15  Format 


NM 

I  NS ( 1 )  *  (NTEHP(I)  -1) 

1-1 


If  this  card  entry  "  0,  omit  and  go  to  (d),  otherwise  go  to  (c). 


f 


(c)  715  Format 

8815  |  jjj  1  NS(Nl)  |  HF(H1,  MS (NT)) 


12  |  10  ]  20 


6E10.3  Format 

1  TF(N1 ,  NS(W1 )  nST-'Wi .  NS(N1  j,WF(Nl ,  ~NS(W1)1  f 


715  Format 


8815  1 


MF(1,  J)  12  10  20 


32 


-nr: 


,1  ■ 


6E10.3  Format 


i 


TF(1, 

j,  1) 

-  TF(1,  j. 

MF(1.  j)) 

• 

• 

• 

715  Format 

8815 

|  N2 

NS(N2) 

MF(N2,  NS(N 2)) 

12 

10 

20~1 

6E10.3  Format 

TF(N2,  NS(N2),  1)  -  TF(N2,  NS(N2),  MF(N2,  NS(N2)) 

N1  Is  first  mission  with  NTEMP(Nl)  -  2 
N2  Is  last  mission  with  NTEMP(N2)  »  2 

(d)  15  Format 


715  Format 

17915  1  1  1  1  |  MLI(l)  1  1  1  12  |  20  1 

6E10.3  Format 
SLI(1,  1)  -  SH(1,  MLI(D) 


715  Format 

7915  1  1  1  1  1  MLI (1 )  1  1  1  12  1  20  1 


6E10.3  Format 


Sll(1,  1)  -  SLI(1,  MLI(O) 


715  Format 


7915  1  1  |  NM  1  MLI(IW)  1  1  1  12  j  20  j 


i 

! 

! 

I 


33 


SLI (NM,  1)  -  SLI (NM,  MLI(l)) 


6E10.3  Format 


SLTP(i,  1)  -  SLTP(1,  HLTP(Q) 


715  Format 


6215  1 

NM  MLTP(NM) 

1  12  20 

6E10.3  Format 

SLTP(NM,  1)  -  : 

SLTP(NM,  MLTP(NM)) 

15  Format 

NH 

715  Format 

|  8515  1  | 

1  |  MLTN(l) 

O 

CM 

CM 

6E10.3  Format 

1 - — - — . . — - - - — - - -  ■  ■■  . 

SLTN(1 »  1)  -  SLTN(1,  MLTN(l)) 


715  Format 

18515  | 

1 

l  1  1  MLTN(i)  |  1  1  12  1  20 

6E10.3  Format 

SLTN(1 , 

1) 

-  SLTN(i,  HLTN(I))  | 

• 

• 

• 

715  Format 

1 851 5 

1 

~1  NM  I  MLTN(NM)  fl  1  12  |  20 

6E10.3  Format 

SLTN(NM,  1)  -  SLTN(NH,  MTLN(NM)) 
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15  Format 


N3  is  the  number  of  missions  that  have  a  touch-and-go  landing. 
If  N3  =  0  omit  this  card  and  go  to  (f),  otherwise  go  to  (e). 


(e)  715  Format 


6E10.3  Format 

STGI(N4,  1)  -STGI(N4,  MTGI ( N4 ) ) 


715  Format 

11215  1  1  1  1  1  HTGI(I)  1  1  1  12  1  20~ 

6E10.3  Format 

STGI ( 1 ,  1)  -  STGI(1,  MTGI(i)) 

715  Format 

11215  \  1  1  N5  |  MTGI(N5)  |  1  j  12  |  20 

6E10.3  Format 

$TGI(N5),  1)  -  STGI(N5,  MTGI(N5)) 


715  Format 


11515 

m 

N4 

MTGG(N4) 

— 

1 

12 

20  j 

16E10.3  Format 

STGG(N4.  1)  -  STGG(N4,  MTGG(N4)) 


715  Format 


|  11515 

1 

1 

MTGG(i) 

1 

12 

20  | 

6E10.3  Format 


STGG(1,  1)  -  STGG(i ,  MTGG(i)) 


715  Format 


|  11515 

1 

N5 

MTGG(N5)  | 

1 

j_J2_ 

20  | 

6E10.3  Format 


STGG(N5,  1)  -  STGG(H5, 

,  MTGG(N5)) 

(f)  15  Format 

NM-NHTMS  | 

715  Format 

|  12460  1  1  1 

NS(1)  3  | 

12  |  10 

1215  Format 

NFF(1,  1,1)-  NFF(1 , 

1.  NS(D) 

• 

• 

• 

715  Format 

■KODiW 

mmwEm 

12  1  10  1 
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1215  Format 


NFF(1 ,  j,  1)  -  NFF(1,  j,  NS(J) ) 


715  Format 


15  Format 


715  Format 
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15  Format 


Q5T) 


715  Format 


115  I  1  I  1  I  MTTP(l)  |  1  |  12 


NTTP(1 ,  1)  -  NTTP(1,  MTTP ( 1 ) ) 


115  1 


MTTP(1 )  1  |  12  20 


NTTP(1,  1)  -  NTTP(i ,  MTTP(i)) 


115  |  1  1  NM  1  MTTP(NM)  j  1  |  12  |  20 


NTTP(NM,  1)  -  NTTP(NH,  MTTP(NM)) 


15  Format 


715  Format 


415  l  1  1  1  MTTN(l)  1  l  12  |  20 


I 

I 


1215  Format 


NFP(i,  j,  1) 

-  NFP(i ,  j,  MFP(i ,  j)) 

• 

• 

• 

715  Format 

715  1  NM 

NS(NM)  MFP(NM,  NS(NM)) 

12 

Ljo 

r*i 

1215  Format 

NFP(NM,  NS{NM),  1)  -  NFP(NM,  NS(NM),  MFP(NH,  NS(NM)) 


15  Format 


NM 

£  NS(1 ) 
i=l 


715  Format 


|  3115  |  1 

1215  Format 

NFN(1,  1,  1) 

-  NFN(1,  1,  MFN(1,  1)) 

• 

• 

• 

715  Format 

1 3115  [  1  1 

j  MFN(1 ,  j)  12  10  20 

1215  Format 

NFN(1,  J.  1)  -  NFN(1,  j,  MFN(1 ,  J)) 


[3115  |  NM  1  NS(NM) 


715  Format 

T  MFN(NM,  NS(NM))  1  12  j  10  j  20~1 
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NFN(NM,  NS(NM) ,  1)  -  NFN(NM,  NS(NM) ,  MFN(NM,  NS(NM)) 


715  Format 


7915  |  1  I  1  ML 1(1)  |  1  I  12  I  20 


1215  Format 


NLI (1 ,  1)  -  NLI(1 ,  MLI(1}) 


715  Format 


7915  1 


1  I  12  20 


1215  Format 


NLI ( 1 .  1)  -  NLI  (1,  MU(D) 


7915  1  1  1  NM  1  MU(NH)  |  1  |  12  |  20 


NLI(NM,  1)  -  NLI(NM,  MLI(NM)) 


8215  |  1  |  1  |  MLTP(l)  |  1  |  12  |  20 


NLTP(1 ,  1)  -  NLTP(1,  NLTP(l)) 


8215 


1 


1 


12 


20 


1215  Format 


1  NLTPCi ,  1)  -  NLTPQ,  MLT_P.(iil 


715  Format 


8215 

1 

NM 

MLTP(NM) 

1 

12 

20 

1215  Format 

NLTP(NM,  1)  -  NLTP(NM,  MLTP(NM)) 


15  Format 

NM 

715  Format 

|  8515  | 

1  1 

MLTN(l)  |  1  12 

20 

1215  Format 

!  NLTN(1 , 

1)  -  NLTN(1, 

MLTN(l)) 

1 

715  Format 

I  8515  | 

1  i  MLTN(i)  |  1  |  12  |  20 

1215  Format 

NLTN(1 , 

1)  -  NLTN{1 i  MLTN(i))  1 

• 

• 

• 

1  8515  | 

1  |  NM  1  MLTN(NM)  1  1  1  12  t  20l 

1215  Format 

NLTN(NM,  1)  -  NLTN(NM,  flLTN(NM)) 
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N1  is  first  mission  with  NTEMP(Nl)  -  2 
N2  is  last  mission  with  NTEMP(N2)  ■  2 
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15  Format 


N3  Is  the  number  of  missions  that  have  a  touch-and-go  landing. 
If  N3  ■  0  omit  this  card  and  go  to  (j)  otherwise  go  to  (i). 


715  Format 

11215  I  1  N4  MTGI(N4)  1  [  12  20 


1215  Format 

NTGI(N4,  1)  -  NTGI (N4,  MTGI ( N4) ) 


715  Format 

11215  ill  MTGI { 1 )  1  12  20 


1215  Format 

NTGI (1 ,  1)  -  NTGI (1 ,  MTGI (1 ) ) 


715  Format 

11215  1  1  1  N5  1  MTGI(N5)  1  \  12  \  20~ 


1215  Format 

NTGKN5.  1)  -  MTGKM5.  MTGI(N5)) 


N4  Is  first  mission  with  NTG(N4)  f  0 
N5  Is  last  mission  with  NTG(N5)  f  0 


15  Format 


P51 


'■SS; 


?-'*  SfirtTifS ,E  i- 
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1215  Format 


MFP(i ,  1)  -  MFP(1 ,  NS(1 ) ) 


715  Format 


11845  I  1  I  NM  I  NS{NM)  1  1  12  10 


1215  Format 


MFP(NM,  1)  -  MFP(NM,  NS(NM)) 


15  Format 

NM 

715  Format 

11965  I  1  1  1  |  N$(l)  j  1  |  12  1  10 


1215  Format 


MFN(1 ,  1)  -  MFN(1,  NS(1)) 


715  Format 


11965 

1 

1 

NS(1 ) 

1 

12 

1215  Format 


MFN(1,  1)  -  MFN(i ,  NS(1 )) 


715  Format 


NM  I  NS(NM)  l  1  1  12  l  10 


1215  Format 


MFN(NM,  1)  -  MFN(NH,  NS(NM 


!■»  -  •; 


15  Format 

B 

715  Format 

12085  1 

i  1 

1  1 

1  NS(1) 

1 

12 

1215  Format 


MFG(1,  1)  -  MFG(1.  NS(1)) 


715  Format 


12085  \  1  |  1  |  NS(1)  i  1  |  12  |  10 


1215  Format 


MFG(1,  1)  -  MFG(1,  NS(1 )) 


715  Format 


12085  1  \  m  \  NS(NM)  \  1  \  12  \  10 


1215  Format 


MFG(NM,  1)  -  MFG(NM,  NS(NM)) 


15  Format 


(NTEMP(I)  -  1) 


If  this  entry  "  0  omit  and  go  to  (1)  otherwise  go  to  (k) 

715  Format 


12250  1  Nl  NS(N1)  |  1  112  |  10 


1215  Format 


MF(N1 ,  1)  -  MF(N1 ,  NS(N1)) 


715  Format 


l 

i 

NS(1 ) 

1 

12 

nri 

1215  Format 

MF(1,  1)  -  MF(1,  NS ( 1 ) )  “ 


715  Format 


[  12250  |  1  \  N2  |  NS(N27 


HJL 


10 


1215  Format 

HF(N2,  1)  -  MF(N2,  NS(N2))  1 

N1  is  first  mission  with  NTEMP(N1  -  2 
N2  is  last  mission  with  NTEMP(N2)  -  2 


(1)  END  OF  FILE 

The  first  card  contains  fourteen  (14)  fixed  point  (Integer) 
numbers  arranged  In  15  fields.  These  entries  are 

(1)  IDENT  -  run  number 

(2)  1 

(3)  1 

(4)  11  ♦  2N  ♦  2L 

N  -  1  if  1  Is  In  [1,  M]  and  NTG(i)  f  0 
N  ■  0  otherwise 
L  ■  1  if  NRANM  -  2 
L  -  0  if  NRANM  -  1 

(5)  0 
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(6)  8  +  M  +  2N 

M  ■  1  If  1  Is  In  [1,  NM]  and  NTEMP(I)  «  2 
H  *  0  otherwise 

(7)  13  +  2M  ♦  2N 

(8)  MONTH 

(9)  DAY 

(10)  YEAR 

(11)  NCASE 

(12)  NM 

(13)  NRANM 

(14)  0 

The  second  card  contains  515  fields  with  the  following  entries 

(1)  NMTMS 

(2)  NFLT 

(3)  NHRS 

(4)  MISRPT 

(5)  MNUMS 

The  third  and  fourth  cards  contain  run  descriptive  Information 
as  desired  by  the  user. 

The  fifth  and  subsequent  cards  are  arranged  as  shown  above. 


5  SAMPLE  PROBLEM 

A  sample  run  Is  Included  to  acquaint  the  user  with  the  Input  deck 
arrangement  and  the  output  listing.  The  case  considered  is  an  aircraft 
with  a  lifetime  of  twelve  hours  that  Is  comprised  of  six  flights.  There 
are  two  missions  and  the  mission  time  and  the  mission  selection  Is  ran- 
dom.  A  temperature  profile  Is  Included  In  mission  1  and  both  missions 
contain  touch-and-go  landings.  The  Input  data  and  the  output  listing 
are  given  below. 
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cr  ritn 


6»TC  t»/I7T» 


rrrrn 


su> 


£l  K*.T  LOW  »»tC  l»u» 


CUriRtNtEft 


unto*  «vrat* 


;»WU> ht  lo»b  and  im>E»>Tint  «t autNct 


SSiR'Rass-fSiS 


f  W  **  •*  *V  <V  n* 


1% - M««OQ 


As  noted  In  the  Introduction,  the  spectrum  printout  as  given  by  this 
routine  may  not  be  useful  for  subsequent  use  In  fatigue  or  fracture 
mechanics  calculation.  It  should  not  be  difficult,  however,  to  modify 
this  program  to  form  the  desired  linkage  through  magnetic  tape  or  other 
device. 
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IV  REFERENCES 


APPENDIX 


SPECF  PROGRAM  LISTING 


This  listing  given  below  Is  a  FORTRAN  extended  language  routine. 
The  DIMENSION  statements  are  set  up  for  an  aircraft  flying  twelve 
missions  with  ten  mission  Intervals  each.  The  EQUIVALENCE  state¬ 
ments  are  compatible  with  a  spectrum  for  fifteen  missions  with 
eight  mission  Intervals  In  each  mission.  Consequently,  this  case 
may  be  analyzed  with  a  simple  DIMENSION  statement  change. 


_ i»a». 


;c i««T35 -  nwaiirw  sris 


4 JO _ rpPW4t  t/10»»  ?QH»tI6HT  J5INJHUH  LQtp  »nCT*U» _ LNADtli 
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